Diving into the proteomic atlas of SARS-CoV-2 infected cells

The COVID-19 pandemic was initiated by the rapid spread of a SARS-CoV-2 strain. Though mainly classified as a respiratory disease, SARS-CoV-2 infects multiple tissues throughout the human body, leading to a wide range of symptoms in patients. To better understand how SARS-CoV-2 affects the proteome from cells with different ontologies, this work generated an infectome atlas of 9 cell models, including cells from brain, blood, digestive system, and adipocyte tissue. Our data shows that SARS-CoV-2 infection mainly trigger dysregulations on proteins related to cellular structure and energy metabolism. Despite these pivotal processes, heterogeneity of infection was also observed, highlighting many proteins and pathways uniquely dysregulated in one cell type or ontological group. These data have been made searchable online via a tool that will permit future submissions of proteomic data (https://reisdeoliveira.shinyapps.io/Infectome_App/) to enrich and expand this knowledgebase.


Results
For this study, all cells analysed were infected "in vitro" with the SARS-CoV-2 virus with the same Multiplicity of Infection (MOI) promoting data homogeneity.Except for SH-SHY5 and HepG2 cells, all other cell lines used in this investigation originated from humans and underwent subsequent differentiation processes, as outlined in the methods section.Before any subsequent analyses, a viral kinetics analysis at three time points indicated that all cells analysed here were indeed infected by SARS-CoV-2.Data regarding T-cells, monocytes, adipocytes, and astrocytes were included in their originally published articles [11][12][13]28 . Desite being able to replicate within CNS cells, the virus was not able to secrete viral particles into the medium; all other cell types showed signs of both viral replication and secretion (Supplementary Fig. 1).
Proteomic analyses of SARS-CoV-2-infected cells were used to build what we will hereby refer to as the SARS-CoV-2 infectome, a compilation of differentially expressed proteins in cells in response to SARS-CoV-2 infection (Fig. 1A, and Supplementary Table 1).Among all analyses, a total of 3098 proteins were quantified, of which 1652 were found differentially regulated (p-value < 0•05, Fig. 1B).Our proteome analysis revealed a high level of data reproducibility and a wide range of identified protein abundance, as corroborated by dynamic range analysis (Supplementary Fig. 2 and 3).It is also noteworthy to emphasize that very few missing values were identified, and less than 2% of the peptides exhibited more than 2 missed cleavages.
Pair-wise Pearson's correlation analysis, followed by hierarchical clustering, showed that almost all cell types clustered as expected by their ontology at the whole-proteome level, a pattern that was markedly distinct while considering only differentially regulated proteins (Supplementary Fig. 4).These differences may be explained by high diversity effects elicited by SARS-CoV-2 in each cell type, which can vary in the number of differentially regulated proteins as well as their regulation direction (Fig. 1C).Interestingly, neuronal cell models were unique in that they presented more upregulated proteins in their infectome.Additionally, we observed that 551 proteins were dysregulated in more than one cell type; however, none were found differentially regulated in all cell types (Fig. 1D, Supplementary Fig. 5).
A pathway enrichment analysis for each infectome identified 151 biological pathways associated with SARS-CoV-2 infection, with the neuronal infectome being the largest contributor (Supplementary Fig. 6).When filtering for pathways enriched in at least 50% of infectomes, SARS-CoV-2 modulated the protein related to energy metabolism (mainly glycolysis), infectious diseases, protein metabolism (both synthesis and degradation processes), signalling/homeostatic pathways, and potential comorbidities (Fig. 1E) in nearly all cells, in line with expected responses to viral infection and viral particle production.Six cellular infectomes (neurons, SH5Y, monocytes, hepatocytes, and both types of adipocytes) also enriched for Coronavirus disease, highlighting the canonical infection pathways in several cell types outside of the respiratory system.A protein-protein interaction network enriched for coronavirus disease reveals several key proteins ubiquitously present dysregulated across all the six cell types (Supplementary Fig. 7A).Analysis of the proteins associated with Coronavirus Disease reveals a deficiency of these pivotal proteins in Astrocytes, Caco-2, and T-Cells, as depicted in Supplementary Fig. 7B.
Given the variety of cellular functions and complexity of the respective infectomes of each cell type, we categorised differential protein expression by ontology: CNS cells (neurons, differentiated SH-SY5Y neuronlike cells, and astrocytes), gastroenterological cells (CACO-2, hepatocytes), white blood cells (monocytes and T-cells), and adipose tissue (differentiated visceral and subcutaneous adipocytes).This allowed an evaluation of similarity among cellular infectomes within ontological groups, useful in increasing our understanding about changes triggered by SARS-CoV-2 to better predict body response to infection.
All results and findings derived from these proteomic analyses were compiled and published in a unique database, forming the SARS-CoV-2 Infectome Atlas.This atlas, being the first of its kind, details the proteomic dysregulations elicited by SARS-CoV-2 in nine cellular infectomes, groupable by their tissue ontology.Individual proteins are searchable, providing information about normalised fold change differences, as well as altered pathways and their constituent proteins.This allows for targeted searching of proteins and pathways of interest within a SARS-CoV-2 context (https:// reisd eoliv eira.shiny apps.io/ Infec tome_ App/).

SARS-CoV-2 infection modulates CNS cell proteomes
The infectome of CNS cells was composed of 810 differentially regulated proteins, with astrocytes presenting the fewest dysregulations (Fig. 2A).Nine differentially regulated proteins were observed in all cell types, primarily associated with translation machinery (RPL29, RPL22, PABPC1, RPL12) and energy metabolism (TKT, PHGDH).Interestingly, however, these proteins presented a distinct regulation profile in each CNS cell, despite similarities in function and role.
In energy metabolism, for instance, PHGDH, a pivotal protein throughout neurodevelopment, was upregulated in all infected cell types; whereas TKT was upregulated in both neuronal cell types but downregulated in astrocytes (Fig. 2B).Astrocytic marker vimentin (VIM), also present in excitatory neurons, was downregulated in infected astrocytes, and upregulated in both infected neuronal cell types (Fig. 2B).This divergent regulation pattern indicates that, even when cells share proximal ontological characteristics, SARS-CoV-2 can induce different and even divergent protein-level responses.
To understand how SARS-CoV-2 can potentially affect these in vitro CNS models at the functional level, the Kyoto Encyclopedia for Genes and Genomes (KEGG) was used for in silico enrichment analyses 32 .Performing a pair-wise comparison of the top ten pathways enriched for each cellular model (Fig. 2C), the NPC-derived neuron infectome was the only CNS cell type seen to incorporate all of the top ten terms found in the other CNS cell types.When both neuronal infectomes were compared, differentiated SH-SY5Y cells enriched 67% fewer terms than NPC-derived neurons, highlighting that SARS-CoV-2 response can differ by cell model, even of the same cell type (Supplementary Fig. 6).Nevertheless, at the pathway level, these cellular infectomes still shared several enriched pathways, such as glycolysis, pentose phosphate pathway, and proteasomes (Fig. 2C).
Since glycolysis was found enriched among the top ten enriched terms of all CNS infectomes (p-adj.< 1•0e-6), we highlighted the proteins differentially regulated in this pathway in each cell type (Fig. 2D).Most of the proteins

SARS-CoV-2 infection modulates gastroenterological cell proteomes
Proteome changes were evaluated in intestinal epithelial cells (CACO-2 cell line) and hepatoma cells (HepG2 cell line), both of which derive from cancer patients and present epithelial-like morphology.Infected CACO-2 cells presented 332 deregulated proteins while HepG2 cells presented 343 deregulated proteins; 38 proteins were deregulated in both cell lines (Fig. 3A).These overlapping proteins are related to vesicle-mediated transport (import into the nucleus, endocytosis, and exocytosis) and cellular metabolic processes (Fig. 3B).Of the terms obtained from KEGG enrichment, 11 pathways were shared between both gastroenterological cell lines (Fig. 3C).Among the enriched terms, we found bacterial invasion of epithelial cells, Salmonella infection, and pathogenic E. coli infection.These processes could be associated with the diarrhoea induced by SARS-CoV-2 infection through a mechanism similar to previously documented bacterial invasion processes in the gut.Among the proteins associated with these pathways, only CDC42, ANXA2 and ARF1 were differentially regulated in both cell types; nevertheless, protein families (such as tubulin or MAPK), are shared in both cell types (Fig. 3D), which could lead to similar effects on epithelial homeostasis.

SARS-CoV-2 infection modulates immune cell proteomes
The immunological cells analysed in our current datasets consist of SARS-CoV-2-infected primary monocytes and T lymphocytes.By testing these two cell types, both myeloid and lymphoid lineages can be tested, covering the innate immune response with monocytes and the adaptive response with the lymphocytes.Our quantitative analysis found 139 proteins dysregulated in T-cells and 463 in monocytes, with 35 proteins dysregulated in both cell types (Fig. 4A).These proteins are associated with cellular metabolic processes; though the heatmap shows that fold change-and therefore likely the cellular response as well-can vary between myeloid and lymphoid cells (Fig. 4B).
Pathways enriched by the monocyte and lymphocyte infectomes (Fig. 4C) presented a large overlap between cell types, of which several terms were associated with immune system responses (such as infections and diseases).Among these terms, we highlight FcγR-mediated phagocytosis, a signalling pathway that interplays adaptive and innate immune responses mediated by antibodies.We also noticed that, within this pathway, lymphocytes presented only upregulated proteins, while monocytes were more diverse in their response (Fig. 4D).Only the monocyte infectome was enriched for COVID-19 disease, suggesting that these myeloid cells may play a pivotal role in the canonically activated pathways during the course of SARS-CoV-2 infection.

SARS-CoV-2 infection modulates visceral and subcutaneous adipocyte proteomes
Since visceral and subcutaneous adipocytes are known targets for SARS-CoV-2 infection and storage 33,34 we investigated the proteomic modulations that occur in these cell types.242 (visceral) and 175 (subcutaneous) proteins were found to be deregulated after SARS-CoV-2 infection.While nearly 60% of identified proteins were www.nature.com/scientificreports/observed in both cell types, 44 proteins were dysregulated in both (Fig. 5A), most of which were in the same direction.No pathways were found to be enriched, despite the presence of many glucose metabolism proteins (Fig. 5B).
When comparing the top enriched pathways from both adipocyte types, only proteasome processes were unique to visceral adipocytes; subcutaneous adipocytes presented more unique processes, namely protein processing in the endoplasmic reticulum, central carbon metabolism in cancer, glucagon signalling pathway, and necroptosis (Fig. 5C).Despite these differences, TCA cycle proteins were found downregulated in both types, indicating dysregulated NAD + and FAD reduction, culminating in lowered ATP production (Fig. 5D).

Discussion
COVID-19 was initially classified as a respiratory disease; however, many studies have since proven that SARS-CoV-2 has the capacity to infect several other cell types and tissues.It has since been classified as a systemic disease with effects including kidney failure, permanent liver damage, diarrhoea, and neurological symptoms and sequelae.An increasing necessity to better understand the molecular and biochemical mechanisms triggered by SARS-CoV-2 infection throughout the human body led us to create the SARS-CoV-2 Infectome Atlas, in which several human cell types infected by SARS-CoV-2 were analysed.www.nature.com/scientificreports/SARS-CoV-2 has been found in human brain tissue collected postmortem and has also been shown to infect CNS models such as organoids and hamster brains 11,35 .Nonetheless, there is still no consensus as to if the brain alterations are directly related to the viral infection or if it is a secondary effect.Moreover, the mechanisms utilised by the virus to infect these cells and which changes are triggered or provoked by its replication are still a matter of discussion.This work reinforces the ability of SARS-CoV-2 to infect different CNS cell types and subsequently trigger proteomic changes.NPC-derived neurons exhibited the most affected model by the size of its infectome.Around 25% of the total identified proteomes from NSC-derived astrocytes, neurons, and differentiated SH-SY5Y neuroblastoma cells were deregulated after infection.Neurons or astrocytes have been hypothesised to be some of the first affected cells within the CNS due to their proximity to the blood-brain barrier (BBB), which would act as the port of entry of the virus into the CNS.Another noncompeting hypothesis is that axonal projections that reach the olfactory bulb may carry the virus into the CNS, thereby bypassing the BBB 5,36 .
Furthermore, infection of the CNS increases IL-6 levels, primarily as a response by astrocytes and microglia, potentially reducing the integrity of the BBB and promoting further viral invasion into the CNS 37,38 .This increase in IL-6 levels has been observed in the cerebrospinal fluid of patients with COVID-19 39,40 .In this study, the astrocyte infectome enriched for biological pathways involved in other chronic inflammatory diseases and responses to infectious agents, suggesting that these cells may be involved in the primary response to SARS-CoV-2 infection.
Within brain tissue, SARS-CoV-2 presence is concentrated in astrocytes despite a more accentuated proteomic change in neurons 11 .While comparing results from the enrichment analysis, NPC-derived neurons and differentiated SH-5YSY cells had a higher overlap of enriched pathways than with astrocytes, likely due to ontological similarity.Many of the enriched pathways were associated with neurodegenerative disease, glucose energy metabolism, COVID-19 disease, and ribosomal translational functions.The enrichment for neurodegenerative diseases could potentially be explained by increases in cell death and apoptosis.
Energy metabolism was affected in all CNS cells analysed in this study, making it a pathway of interest for further study.Proteins belonging to glycolytic pathways were generally upregulated in neurons and SH-5YSY cells, while in astrocytes they were downregulated.Calmodulin-dependent protein kinase (CAMK2D) was found upregulated in astrocytes, a protein that plays an important role in glutamatergic synapses through the post NMDA receptor activation events, possibly representing a higher uptake of glutamate from the extracellular medium.
In all cell types, PHGDH was found upregulated (Fig. 2B), a protein which catalyses the synthesis of L-serine, the precursor to D-serine, which is essential for NMDA receptor function throughout neurodevelopment.The upregulation of PHGDH in these cells could be related to an increased consumption of glutamate by the infected CNS cells to offset the increased energy demands resulting from the replication of viral particles.Crunfli et al. showed an imbalance in the glutamate/glutamine interplay between SARS-CoV-2 infected astrocytes and neurons 11 and proposed that the higher mortality observed in neurons and SH-5YSY cells after being exposed to medium conditioned with infected astrocytes is due to a disruption of this feedback system between astrocytes and neurons.
More than half of COVID-19 patients report gastrointestinal symptoms, most commonly diarrhoea, but also including anorexia, nausea, vomiting, abdominal pain and gastrointestinal bleeding 41 .ACE2 is abundantly expressed in the intestinal epithelial cell cytoplasm 42 , enabling SARS-CoV-2 to directly target gastrointestinal cells, specifically gastric and intestinal epithelial cells 43 .Based on this evidence and the various gastroenterological symptoms, we chose to investigate the cellular response to infection in epithelial intestinal cells (CACO-2) and hepatocytes when creating our SARS-CoV-2 infection proteome map.The pathway analysis from the 332 proteins in the CACO-2 infectome mainly enriched for pathways related to chronic inflammatory diseases, energy metabolism, and biosynthesis of amino acids, highlighting the necessity of SARS-CoV-2 to induce changes in energy and protein metabolism to quickly replicate and create new viral particles.
Bojkova et al. used the same cell line labelled with stable isotopic amino acids (SILAC) to study the proteomic profile after SARS-CoV-2 infection 44 , discovering 459 proteins dysregulated by viral infection, enriching for glycolysis and nucleotide metabolism pathways.They also tested inhibitors of protein translation, RNA splicing, glycolysis, and nucleotide synthesis at non-toxic concentrations, which, in line with our findings, inhibited viral replication.Another study raised the hypothesis that, upon the ACE2-mediated entry of SARS-CoV-2 into the gastrointestinal tract via the small intestine, viral invasion and expansion occur, triggering gastrointestinal inflammation and potentially evolving into a strong cytokine response 45 .Within the gastroenterological system, the liver is the organ that is considered to be most impacted by COVID-19, with approximately 58-75% of patients presenting hepatic injury and/or increased levels of transaminases [46][47][48][49] .While patients with COVID-19 generally present only moderate liver damage, more severe cases have been reported 47,48 .Wang and colleagues were the first to provide evidence of cytopathy in hepatocytes caused by SARS-CoV-2 infection that could cause liver impairment in the infected host 49 .Our analysis indeed shows that SARS-CoV-2 elicits proteomic changes in the HepG2 cell line, affecting proteins related to energy metabolism, the proteasome, amino acid biosynthesis, and pyruvate metabolism.When filtering for proteins associated with carbon metabolism in hepatocytes, there was an increase in proteins involved in pyruvate production, indicating once again that glucose metabolism is crucial during viral replication.
Adipose tissue has also been associated with SARS-CoV-2 infection and the subsequent inflammatory response 34,[50][51][52] .Since the entry of SARS-CoV-2 into the digestive system, obesity, and high quantities of adipose tissue are risk factors and indicators of poor prognosis for severe COVID-19 and its evolution 53 , we evaluated infectomes of primary stromal-vascular cells isolated from subcutaneous (abdomen) or visceral (omentum) adipose tissues.Subcutaneous and visceral adipose tissue cells presented different responses to infection, with only 44 proteins in common.A pathway enrichment analysis highlighted changes in pyruvate metabolism, carbon metabolism, glycolysis/gluconeogenesis, and the TCA cycle in both cell types.www.nature.com/scientificreports/These dysregulations could be associated with the hyperglycemia that is often reported in acute COVID-19 patients and other metabolic dysfunctions that is seen in this disease 54 .Hospitalised COVID-19 patients with hyperglycemia also presented a higher prevalence of insulin resistance when compared to ARDS (acute respiratory distress syndrome)-positive and ARDS-negative controls, due to changes in adipokine levels (adiponectin and leptin), instead of beta-cell function 54 .This was also observed in other studies that analysed adipokine levels in infected hamsters 55,56 .The energy metabolism dysfunctions seen in both adipocyte infectomes could be disrupting the proper functioning of these cells, inducing changes in adipokine secretion.When epithelial intestinal cells (CACO-2), hepatocytes (HepG2), and both subcutaneous and visceral adipose tissue were compared, only MDH2 (malate dehydrogenase 2) was found in all infectomes, being upregulated in CACO-2 cells and downregulated in all other cell types (Supplementary Fig. 8).
The immune response that is mounted upon SARS-CoV-2 infection is against the spike and nucleocapsid viral components and seeks to eliminate the virus from the host 57 .Lower levels of neutralising antibodies in circulation and a progression to severe COVID-19 are strongly associated with an exacerbated inflammatory response, alterations in the immune system response characterised by protein synthesis reduction, a "cytokine storm", lymphocytopenia, and T-cell exhaustion 58,59 .Immune responses occur through two different mechanisms: the innate immune response, mediated by macrophages, dendritic cells, and natural killer cells, and/or the adaptive immune response, mediated by T-cells and B-cells.
To cover as many of the proteomic changes brought about by SARS-CoV-2 infection as possible, we used human primary T-cells and monocytes as infection models, thereby screening both innate and adaptive responses 8,12,28 .The resulting pathway analyses showed a high similarity between the T-cells and monocytes, enriching for several pathways related to response to infection, phagocytosis, neutrophil cellular trap formation, and, once again, energy metabolism.Ribosomal changes and coronavirus disease were identified only in monocytes, likely since monocytes are responsible for the activation of macrophages and produce the initial immune response.In agreement with this hypothesis, chemokine signalling activation and leukocyte transendothelial migration were both enriched exclusively in monocytes.The T-cell infectome also exclusively enriched for certain biological pathways, which included apoptotic and necroptotic pathways, supporting a study that saw a decrease in CD4 + T-cell viability after SARS-CoV-2 infection 12 .As SARS-CoV-2 replication is dependent on the presence of glucose, this highlights the many pathways involved in glucose metabolism that were enriched in both cell types.Moreover, this suggests that energy production dysfunctions could be directly related to an exacerbated immune response and resulting cytokine storm 28 .
When all cell types were compared together, a large overlap was observed in identified proteins; however, when only dysregulated proteins were considered, this overlap vanished (Supplementary Fig. 5).Nie et al. observed a similar trend in their proteomic analysis of seven tissue types from patients who died from COVID-19 8 .They showed that only 39 proteins were commonly dysregulated of a total of 11,394 identified proteins, representing less than 0•004% of identifications.Despite the lack of overlap of individual proteins, energy metabolism was affected in every infectome (Supplementary Fig. 9), which strongly supports the hypothesis that glucose metabolism is a key factor for SARS-CoV-2 replication inside the host cell.This hypothesis is further supported by several studies that have reported disturbances in energy production in tissues and cells infected by SARS-CoV-2 8,11,28 .
To facilitate the future investigation of all the data collected and analysed in this work, we have constructed an open-source tool to facilitate the search for proteins found to be dysregulated in each cell type and tissue type.Other data include their relative abundance in comparison to controls and the biological pathways they enrich, all of which is freely available at https:// reisd eoliv eira.shiny apps.io/ Infec tome_ App/.Future proteomic analyses can then be submitted to this database to expand this SARS-CoV-2 Infectome Atlas.
All products used for cell culture are from Thermo Fisher Scientific, MA, USA.The characterization of the BR-1 lineage as astrocyte cells has been previously described elsewhere [14][15][16][17] .We generated eight batches of human astrocytes from BR-1-derived NSCs.The NSCs were of different passages and were used as biological replicates in independent experiments.Our internal control showed that approximately 97% of the neural stem cell-derived astrocytes in culture expressed GFAP, 80•4% expressed vimentin, and 12•9% expressed SOX-2 (markers of progenitor cells).The neural stem cell-derived astrocyte culture expressed more astrocytic markers than Vol:.(1234567890 www.nature.com/scientificreports/progenitor cell markers, indicating an excellent effectiveness of human astrocyte generation.Method described in our previous work 11 .
NSCs were plated on Geltrex-coated plates and maintained with NEM medium at 37 °C in humidified air with 5% atmospheric CO 2 .Upon reaching 40% confluency, the medium was changed to neuronal differentiation medium consisting of 1:1 DMEM/F12: neurobasal medium with 1% B27 supplement (Thermo Fisher Scientific, Carlsbad, CA, USA) and 1% GlutaMAX (Thermo Fisher Scientific, Carlsbad, CA, USA).The medium was renewed every four days over the course of 20 days by removing half of the volume and replacing the same volume of fresh medium.Medium renewal was performed in this manner since factors secreted by the differentiating cells are important for successful differentiation.Two control cell lines were used: GM23279A, obtained from a female subject (available at Coriell; RRID:CVCL_F178) and BR-1 (RRID:CVCL_C062) 14 .Both cell lines were cultivated following previously published protocols 16,[24][25][26] and have been extensively characterised elsewhere [21][22][23] .We also used FACS to analyse cellular markers and found a bona fide neuronal phenotype with the expression of the neuronal markers synaptophysin (75•9% of cells), MAP2 (99•9%), and β-tubulin (99%), as well as astrocytic marker GFAP (8•1%).

Human adipose tissue mesenchymal stem cell isolation, culture and adipocyte differentiation
Patients, samples collection and adipocyte differentiation were previous described in our work 13 .For this work it was used the same samples but with independent runs in the LC-MS/MS system.Briefly, human adipose tissue-derived mesenchymal stem cells (hADSCs) were isolated from abdominal subcutaneous adipose tissue and visceral omental adipose tissue of three individuals who underwent abdominal surgery (i.e., bariatric surgery or cholecystectomy) at the Clinical Hospital of the University of Campinas, Campinas, SP, Brazil.These cells were isolated prior to the COVID-19 pandemic and hence the donors were not infected with SARS-CoV-2.The subjects received written and oral information before providing written informed consent for the collection of the biopsy and use of the tissues.All material used in the procedure was sterile.Biopsies were collected during surgery and transported to the laboratory within 30 min in sealed, sterile falcon tubes for initiation of the procedure of stromal-vascular fraction isolation.The tissue was weighed, cut into small pieces on a petri dish, and digested with 25-30 mL of lysis buffer (1 mg/mL collagenase in HBSS containing 2% BSA, filtered through a 0•22 μm filter) at 37 °C for 30-50 min with slight agitation, until homogeneous.The homogenate was filtered through a 250 μm filter and collected in a sterile falcon tube.After a brief rest period (approximately 5 min), the infranatant containing the stromal-vascular fraction was collected using a Pasteur pipette and centrifuged for 5 min at 200×g and 4 °C.The supernatant was discarded, and the cell pellet was washed with HBSS.The centrifugation and washing steps were repeated twice.The pellet was then resuspended in BM-1 medium (ZenBio) with 10% foetal bovine serum (FBS) and 1% penicillin/streptomycin (penicillin 10,000 units/mL, streptomycin 10 mg/mL).
Cells were cultured at 37 °C and 5% atmospheric CO 2 until semi-confluent (80-90%).They were then trypsinized and frozen in freezing medium (10% DMSO, 50% FBS, and 40% basal culture medium) and stored in a liquid nitrogen biorepository for further assays (e.g., adipocyte differentiation followed by exposure with agents, lipid quantification, determination of cell viability, gene expression analysis, proteomics, immunodetection of proteins, quantification of cell secretion products, and metabolic analysis).

Virus strain
The HIAE-02-SARS-CoV-2/SP02/human/2020/BRA (GenBank accession number MT126808•1) virus strain was used for all in vitro experiments.The virus was isolated from the first confirmed case of COVID-19 in Brazil and kindly donated by Prof. Dr. Edison Durigon (ICB-USP).VSV-eGFP-SARS-CoV-2 was engineered and donated by Prof. Dr. Sean P.J. Whelan (Department of Medicine, Washington University School of Medicine, St. Louis, MO, USA) for SARS-CoV-2 entry experiments 27 .Viral stock was propagated in Vero CCL-81 cells (ATCC,

Blood sample collection and lymphocyte separation
Each COVID-19 patient had heparin and plain blood tubes collected.Whole blood, serum, and plasma samples were separated.Peripheral blood mononuclear cells (PBMCs) from patients and buffy coats were obtained via Histopaque-1077 density gradient (Sigma-Aldrich).Samples were diluted in Hanks balanced salt solution (1:1) and gently poured into 15-or 50-mL conical tubes containing 3 or 10 mL of Histopaque, respectively.Samples were centrifuged at 1400 rpm for 30 min at 4 °C without acceleration or braking.After the PBMC layer was transferred to new tube, lymphocytes were sorted and incubated overnight with RPMI 1640 (Gibco) containing 10% foetal bovine serum (FBS) and 1% penicillin-streptomycin (P/S), at 37 °C with 5% atmospheric CO 2 .Is this work it was reanalalyzed the proteomic data obtained by our previous work.Patients steatments and ethical issues are also described at Davanzo et al. 12 .

In vitro astrocyte infection
Astrocytes were infected with SARS-CoV-2 for 1 h (MOI of 0•1 for proteomic, viral load, and bioenergetic assays) under gentle agitation at room temperature.After viral adsorption, cells were washed twice with phosphatebuffered saline (PBS) and incubated with DMEM/F12 supplemented with 10% FBS, 1% GlutaMAX, and 1% penicillin and streptomycin for 24 h at 37ºC and 5% atmospheric CO 2 .

LC-MS/MS sample preparation, analyses, and data processing
All cellular experiments underwent the same preparation process, as follows.Cells were lysed chemically with lysis buffer (100 mM Tris-HCl, 1 mM EDTA, 150 mM NaCl, 1% Triton-X, and protease and phosphatase inhibitors) then mechanically with an ultrasonication probe during three cycles of 20 s with 90% frequency, on ice.The total protein extract was quantified by Pierce BCA protein assay (Thermo Fisher Scientific, MA, USA) according to the manufacturer's instructions.30 µg of total protein extract from each sample was transferred to a Microcon-10 centrifugal filter with a 10 kDa cutoff for FASP protein digestion 29 .In brief, proteins were reduced (10 mM dithiothreitol) at 56 °C for 40 min, alkylated (50 mM iodoacetamide) at room temperature for 20 min in the dark, and digested overnight with trypsin at 37 °C in 50 mM ammonium bicarbonate (AmBic), pH 8•0.Peptides were recovered from the filter in 50 mM AmBic and trypsin activity was quenched by adding formic acid (FA) to a final concentration of 1% (v/v), whereupon the peptides were desiccated in a SpeedVac vacuum concentrator and stored at -80 °C until analysis.
Digested peptides were resuspended in 0•1% FA.LC-MS/MS analyses were performed in an ACQUITY UPLC M-Class System (Waters Corporation, Milford, MA) coupled online to a Synapt G2-Si mass spectrometer (Waters Corporation, Milford, MA). 1 μg of peptides were loaded onto a trapping column (Symmetry C18 5 μm, 180 μm × 20 mm, Waters Corporation, Milford, MA) and subsequently separated in the analytical column (HSS T3 C18 1•8 μm, 75 μm × 150 mm; Waters Corporation, Milford, MA).For gradient elutions, 0•1% FA was used as eluent A and acetonitrile-FA (99•9% ACN:0•1% FA) as eluent B. A reversed phase gradient was carried out over a 120-min method, with a linear gradient running from 3 to 60% eluent B over 90 min at 300 nL/min.In the Synapt G2-Si, the peptide spectra were acquired by ion mobility-enhanced data-independent acquisition (HDMS E ).Mass spectrometry analysis was performed in Resolution Mode, switching between low (4 eV) and high (25-60 eV) collision energies, using a scan time of The raw data from each experiment were processed in Progenesis QI for Proteomics, version 4•0•× (Waters Corporation, Milford, MA).Tandem mass spectra were searched against a reviewed Homo sapiens proteome (UNIPROT, release 2020-04), using tolerance parameters of 20 ppm for precursor ions and 10 ppm for product ions.For peptide identification, carbamidomethylation of cysteine was set as a fixed modification and oxidation of methionine as a variable modification, up to two missed cleavages were permitted and the false discovery rate (FDR) was limited to 1%.Protein identification was performed using a minimum of one fragment ion per peptide, three fragment ions per protein, and one peptide per protein.Sample acquisition until this point was performed by experimenter blind to group assignment and outcome assessment.Analyses after this point were performed with predefined algorithms and cutoffs, reducing the possibility of experimenter-induced bias.
The label-free quantitative analysis was carried out using the relative abundance intensity of all peptides of a protein with at least one unique peptide after normalisation by all peptide intensities.The expression analysis was performed considering the technical replicates for each experimental condition, following the hypothesis that each group is independent.Proteins with ANOVA (p) ≤ 0•05 between groups were considered differentially expressed.Data integration and visualization was carried out in OmicScope package 62 (v.1.2.2), using Progenesis Method and default parameters.The method applied for proteomics propose are better described at Crunfli et al. 11 .

RNA extraction and viral load
Total RNA extraction was performed using TRI Reagent (Sigma, St Louis, USA) according to the manufacturer's instructions.RNA concentration was determined by a DeNovix spectrophotometer and RNA integrity was assessed by visualisation of 28S and 18S ribosomal RNA on a 1% agarose gel.Reverse transcription was performed with 0•5 µg of RNA using a GoScript reverse transcriptase kit (Promega, Madison, WI, USA) according to the manufacturer's instructions.qPCR was performed using astrocyte cDNA diluted 1:10 and the qPCR SybrGreen Supermix (Qiagen, Valencia, CA, USA) containing forward and reverse primers in RNAse-free water.All reactions were performed in a CFX384 Touch real-time PCR detection system (Biorad, Hercules, CA, USA) and cycling conditions were set as follows: 50 °C for 2 min; 95 °C for 10 min; (95 °C for 15 s; 60 °C for 1 min) × 40 cycles.To evaluate primer specificity, a melting curve analysis was performed by heating samples from 65 to 99 °C (1 °C increment changes at 5 s intervals).All sample measurements were performed in duplicate.Primers were designed with PrimerBlast and used at a concentration of 200 nM.Data were normalised to the expression of 18S (Fwd 5′ CCC AAC TTC TTA GAG GGA CAAG 3′; Rev 5′ CAT CTA AGG GCA TCA CAG ACC 3′) and the relative quantification value of each target gene was determined using a comparative CT method 30 .For virus detection, SARS-CoV-2 nucleocapsid N1 primers (Fwd 5′ CAA TGC TGC AAT CGT GCT AC 3′; Rev 5′ GTT GCG ACT ACG TGA TGA GG 3′) were used as previously described 28,31 .Data were expressed as mean ± SEM.Statistical significance was calculated by a two-tailed unpaired Student's t-test.All analyses were performed using GraphPad Prism 8•0 (San Diego, CA, USA) and a significance level of p ≤ 0•05 was adopted.

Ethical approval
The protocols for human sample collection, storage, and analysis were conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Research Ethics Committees of the University of São Paulo and University of Campinas (CAAE 48,836,721.3.0000.5440and CAAE 78,577,417.8.0000.5404).The latter was confirmed by the National Ethics and Research Council.Procedures for sample collection, processing, and subject information are described below.Donors were not compensated.

Figure 1 .
Figure 1.SARS-CoV-2 infects cells from different tissues and elicits changes in protein abundance.(A) Cell lines from distinct tissues and biological systems were infected with SARS-CoV-2 in vitro and the infectomes (proteins differentially regulated due to viral response) were obtained by shotgun proteomics.(B) Total number of identified proteins (gray) and differentially regulated proteins (colored) in each cell type.(C) Proteins differentially regulated in the respective infectomes separated by direction of regulation.Bubble size indicates the number of differentially expressed proteins (DEPs).(D) Overview of SARS-CoV-2 infectome in all cell types (outer circle), showing proteins that are shared among groups as well as their respective regulation in the heatmap (inner circle).(E) Sankey plot of pathways enriched in at least 50% of infectomes.

Figure 2 .
Figure 2. SARS-CoV-2 induces changes in the proteomes of CNS cell types.(A) Venn diagram showing differentially expressed proteins among NPC-derived astrocytes and neurons and differentiated SH-SY5Y neurons.(B) Normalized abundance of the differentially regulated proteins associated with energy metabolism (PHGDH, TKT) and cellular marker (VIM) for CNS cells.Gray bars represent the abundance in uninfected cells.(C) Pathway enrichment analyses for each CNS infectome.The top 10 pathways were selected for each cellular model, and then a pair-wise comparison was performed to visualize enrichment in other cell types.Bubble size indicates adjusted p-value on a − log10 scale.(D) Upregulated (red) and downregulated (blue) proteins associated with glycolysis in different CNS cell types.Proteins in bold and italic are pivotal in glycolysis and lactate generation, respectively.

Figure 3 .
Figure 3. SARS-CoV-2 induces changes in the proteomes of gastroenterological cell types.(A) Venn diagram showing differentially expressed proteins in HepG2 and CACO-2 cells.(B) Heatmap of fold changes for each protein found in both gastroenterological cell types.(C) Pathway enrichment analyses for hepatocytes and CACO-2 cells.Bubble size indicates adjusted p-value on a − log10 scale.(D) Upregulated (red) and downregulated (blue) proteins related to bacterial invasion of epithelial cells.Proteins highlighted in bold are from the MAPK and tubulin families.

Figure 4 .
Figure 4. SARS-CoV-2 induces changes in the proteomes of immunological cells.(A) Venn diagram showing differentially expressed proteins in monocytes and lymphocytes.(B) Heatmap of fold changes for proteins found in both immune cells.(C) Pair-wise comparison of top pathways enriched against the KEGG database.Bubble size indicates adjusted p-value on a − log10 scale.(D) Upregulated (red) and downregulated (blue) proteins involved in Fc gamma R-mediated phagocytosis.

Figure 5 .
Figure 5. SARS-CoV-2 induces changes in the proteome of adipocytes.(A) Venn diagram showing differentially expressed proteins in visceral and subcutaneous adipocytes.(B) Heatmap of fold changes for proteins found in both adipocyte types.(C) Pair-wise comparison of top pathways enriched against KEGG database.Bubble size indicates adjusted p-value on a − log10 scale.(D) Upregulated (red) and downregulated (blue) proteins involved in the citrate (TCA) cycle. https://doi.org/10.1038/s41598-024-56328-3 for ion mobility separation was 1000 m/s and the transfer wave velocity was 175 m/s.A human [Glu1]-Fibrinopeptide B standard (Waters Corporation, Milford, MA) was used as the reference lock-mass compound.Each sample was acquired in at least duplicate.